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A  TWO-DIMENSIONAL  NUMERICAL  PROGRAM 
SIMULATING  THE  EFFECTS  OF  FREE-ELECTRON  LASER 
NONROTATIONAL  SYMMETRY 

(Romanized  title:  *Ziyou  Dianzi  Jiguang  Feixuanzhuan  Duichen  Xiaoying 

De  Envei  Shuzhi  Monf) 

by  Qian  Sihai  and  Shi  Yijin 

(China  Institute  of  Atomic  Energy,  P.O.  Box  275  (18),  Beijing  102413) 


Abstract:  A  nonrotational  symmetry  2D  FEL  numerical  simulation  program  was  set 
up  with  electron  wiggling  motion  direction  extended  infinitely.  This  program  was  used  to 
calculate  light  field  evolution  under  conditions  of  nonideal  incidence.  Results  showed  that 
both  shifting  the  center  of  the  electron  beam  and  tilting  the  beam  axis  caused  the 
exponential  degradation  of  luminous  power  output  and  excitation  of  higher-order  transverse 
modes,  as  well  as  deviation  of  the  light  beam’s  baiycenter.  In  addition  to  emittance, 
transverse  phase  space  geometry  of  the  electron  beam  is  also  an  important  parameter  of 
FEL  operation  quality  and  has  an  influence  on  luminous  power. 

Key  terms:  nonrotational  symmetry,  nonideal  incidence,  phase  space  geometry 
1.  Introduction 

The  higher  dimeBsion  effect  of  free-electron  lasers  [FELs]  is  a  problem  that 
[scientists]  have  wanted  to  understand  for  a  long  time  but  have  not  yet  fully  understood. 
It  is  not  only  important  that  its  effects  on  transverse  mode  competition  and  light  guide 
phenomena  be  understood  for  the  field  of  physics,  but  the  field  of  engineering  design  should 
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also  consider  its  effects  on  gain,  extraction  efficiency,  and  light  beam  quality.  Because  there 
is  a  larger  amount  of  computation  work  for  the  PEL  3D  [program],  it  is  necessary  to 
simplify  the  model  even  more  and  reduce  the  quantity  of  computations  of  individual  plans 
in  order  to  clearly  understand  the  essence  of  this  problem.  Considering  [this  problem]  as 
a  two-dimensional  approximation  can  decrease  computation  time  by  at  least  one  magnitude. 
Thus,  many  people  have  obtained  good  results  by  using  rotational  symmetry 
approximations^^'^'  and  reducing  the  problem  to  two  dimensions.  Even  though  the 
geometrical  structure  of  FELs  is  ordinarily  close  to  rotational  symmetry  (and  rotational 
symmetry  approximations  are  thus  important),  there  are  still  many  situations  that  give  rise 
to  nonrotational  symmetry.  As  far  as  the  problem  of  coUimation  of  the  three  axes  in  the 
PEL  (the  electron  beam  axis,  the  optical  axis,  and  the  oscillator^  magnetic  axis)  is 
concerned,  it  is  impossible  for  engineering  design  and  physics  debugging  to  completely  avoid 
small  deflections  between  the  axes.  We  can  only  make  demands  for  an  allowable  amount 
of  deviation. 

Because  of  the  special  geometric  structure  of  the  linearly  polarized  oscillator,  wiggler 
field  transverse  distribution  is  unavoidably  elliptical.  [We]  cannot  rule  out  carrying  out 
investigations  that  fully  use  the  elliptical  nature  of  the  wiggler  field  to  inject  oval  cross- 
section  electron  beams.  In  addition,  because  of  the  introduction  of  glancing  incidence  lenses 
into  the  ring-shaped  cavity,  the  light  beam  basically  becomes  an  oval  cross  section  (because 
the  light  waists  of  two  directions  do  not  coincide).  Thus,  asymmetry  exists  in  two  directions 
in  the  transverse  direction.  A  geometrical  approximation  like  this,  that  is,  another  kind  of 
two-dimensional  approximation  relative  to  rotational  symmetry,  is  to  extend  one  direction 
infinitely  so  that  it  becomes  translation-invariant.  Here,  the  direction  chosen  is  the  one 
perpendicular  to  the  wiggler  field,  which  is  called  the  x  direction.  The  wiggler  field 
[direction]  is  the  y  direction,  and  the  electron  longitudinal  motion  direction  is  the  z 
direction. 

In  this  article,  we  observed  and  studied  two  situations  of  nonideal  electron  beam 
incidence;  in  one,  the  electron  beam  has  a  slight  deviation  from  the  oscillator  when  injected 
towards  the  center,  and  in  the  other,  the  electron  beam  incidence  direction  has  a  small  angle 
of  deflection  from  the  oscillator  axis. 


The  word  bozhenqi,  which  I  translated  here  as  "osciUator,"  was  not  found  in  the  dictionary.  A  more  literal 
translation  would  be  "wave  vibrator."  Zhendangqi  is  the  most  common  word  for  oscillator. 
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Of  these  equations,  numbers  (1)  and  (2)  are  electron  longitudinal  motion  equations,  (3)  and 
(4)  are  transverse  motion  equations,  (5)  is  a  complex  light  field  E(y,  z)  =  E(y,  z)exp{/0(y, 
z)}  Maxwell  equation  adaxial  form,  Y  =z/L^  is  the  overall  length  of  the  oscillator,  kj^z) 
=  InfXJ^z)  and  k^^liclX^  are,  respectively,  oscillator  wave  number  and  light  field  wave 
number,  and  ajy,  z)  is  the  reduced  vector-potential  of  the  oscillator. 


a,(y  .2)  = 


me 


k^{z) 


me 


-c.Mli+j  (k.yf] 


^  k.a\{z) 
8y^A:,(z) 


(6) 

(7) 


/q  and  /i  are,  respectively,  zero-  and  first-level  Bessel  functions.  <...  >  indicates  an  average 
of  the  corresponding  amount  of  electrons  in  the  grids  with  difference  grid  coordinates  of  (y, 
z).  Electric  current  intensity  /  in  formula  (5)  is  expressed  in  units  of  A/cm^. 

The  electron  motion  equation  was  derived  through  a  fourth-level  Runge-Kutta 
method,  and  the  light  field  equation  was  derived  through  a  steady-state  difference  format. 

The  initial  values  of  the  electron  motion  equation  (X,  ijr,  y,  and  fi^)  are  usually 
sampled  and  given  m  four-dimensional  phase  space.  Experience  shows  that  the  initial  value 
of  ^  is  very  sensitive  to  calculation,  and  that  there  are  large  discrepancies  between  different 
distribution  sampling  results.  This  is,  of  course,  caused  by  there  being  an  insufficient 
munber  of  sampled  electrons.  Obviously,  increasing  the  number  of  electrons  N  in  the 
sample  in  order  to  stabilize  the  calculation  is  not  a  sound  method.  Because  ijr  brought  about 
serious  consequences  in  fom  variable  quantities,  we  adopted  a  stationary  initiation  method^^. 
We  split  up  a  number  of  samples  N  of  electrons  into  m  groups,  each  containing  an  equal 
number  n  of  electrons.  N  =  mn,  and  the  number  of  electrons  n  in  each  group  all  had  the 
same  (y,y,  and  j5y)  initial  sampling  values,  while  each  group’s  phase  was  fixed  at  ^i=2jri7n, 
i  =  1,  2, ...,  n.  By  choosing  m  =512,  n =8,  we  obtained  stable  results. 

/351 

We  used  two  methods  to  calculate  sampling  efficient^.  One  was  from  evolution  of 
optical  power  r)Q={P{z)-P-^f^YSS^)>  and  the  other  was  from  electron  energy  deficit 
»7e=?[^i(0)-yi(^)]/yi(0)-  During  the  entire  calculation,  the  relationship  TjQ{z)=tjJ^z)  was 
strictly  kept.  Thus,  the  calculation  was  reliable  and  self-consistent. 
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Table  1  Simulation  Parameters 


.  -  .  ^ _ 

electron  bcara 

undulator 

Lorentz  factor 

y  2=  100.32 

Length 

500cm 

Energy  spread 

hh  •/• 

X  =  8cin 

Emittance 

«  =>'^,-4.4  X  10‘*cm -rad 

Period 

Beam  radius  (adjustable) 

*  . — ^0-7854cm  ' 

Current 

/=4.2kA/cm 

Peak  magnetic  field 

B.=  2430Gs 

Optical  field 

Number  of  sampling  electrons  A^*=4096 

Wavelength 

lO.bMin 

Transverse  width  y, 

^  =  3cm 

Input  power 

IMW 

Input  waist  size 

0.35cm 

To  make  a  comparison  with  other  calculations^^'^^  and  to  test  and  verify  the  reliability 
of  the  program,  see  Table  1  for  the  parameters^^^  used  throughout  the  calculation.  Electric 
current  alone  was  found  to  be  4244A,  which  was  times  that  of  document  [2] 

(r= 0.3cm).  The  subjective  consideration  that  prompted  seeking  this  ratio  was  to  cause  the 
electric  current  densities  of  the  two  to  be  equal.  Therefore,  for  the  output  optical  power 
between  the  two,  [their]  power  densities  should  also  be  compared  (2rj7cr^  is  the  ratio  of 
the  rectangular  cross  section  to  the  circular  cross  section,  and  the  X  direction  takes  unit 
length).  Figure  1  shows  output  optical  power  at  different  amounts  of  incident  electron  beam 
center  energy  y.  When  y  =  100,  density  of  output  optical  power  is  IGW/cm^.  Thus,  average 
power  density  is  consistent  with  document  [2]. 

3.  Matching  of  electron  beam  emittance  and  oscillator  output 


Electron  beam  emittance  is  an  important  parameter  of  free  electron  laser  operation 
quality.  The  emittance  of  short  wavelength  free  electron  lasers  especially  requires 
wavelength  magnitude.  Accordingly,  the  electron  beam  should  realize  full  coincidence  with 
the  light  beam  to  the  greatest  extent  possible.  This  was  suggested  through  experience  and 
theoretical  calculations.  From  the  point  of  view  of  electron  beam  transport,  so-called  "full 
coincidence^  means  that  the  electron  beam  envelope  evolves  according  to  the  evolution  of 
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the  light  beam  envelope.  For  electron  beam  transmission,  there  are  not  only  requirements 
for  beam  waist  size  and  position,  the  largest  envelope  radius,  or  other  individual 
measurements,  but  for  the  entire  transmission  process.  Thus,  requirements  must  be  proposed 
not  only  for  emittance,  but  for  the  initial  phase  geometry  of  electron  beams  as  well. 


p«>wcr  W 


Fig.  1  Output  power  at  different  levels  of  incident  electron  beam  energy  y 

For  the  transport  system  of  transverse  motion  determined  in  formulas  (3)  and  (4) 
(including  alternating  quadrupole  focusing),  it  is  easy  to  establish  the  envelope  evolution 
equation  of  the  electron  beam  within  the  system  and  to  solve  it.  The  envelope  radius 
expression  formula  for  ideal  injection  conditions  is 

R\z)  =  R\coiKk,2)+  -^^sin(2V) +  (/?,'*+  ^  (8) 

/352 

Here,  ideal  injection  means  that  the  electron  beam  positioning  and  orientation  have  no 
deviation,  and  that  the  injection  phase  oval  center  is  coincident  with  the  phase  space 
coordinate  center.  (Nonideal  injection  is  discussed  in  the  following  section.)  Qearly,  it  is 
not  only  emittance  e  that  determines  evolution  R{z),  but  also  electron  beam  initial  radius 
/?o  and  slope  Rq  -dRJdz,  where  k^=k^J^/lr  is  the  betatron  motion  wave  number,  which 
is  determined  by  the  intensity  of  the  focusing  field. 
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Considering  the  effects  of  the  electronic  envelope  on  the  evolution  of  the  light  beam 
envelope  and  finding  the  conditions  of  fullest  coincidence  becomes  a  very  complex  nonlinear 
process,  especially  when  gain  is  relatively  large. 

To  emphasize  the  importance  of  incident  phase  space  geometry,  obtain  a  simpler 
situation,  and  obtain  Rq  =0  (the  electron  beam  waist  is  obtained  at  the  oscillator  entrance), 
here,  formula  (8)  can  be  simplified  to 

/t'(2)=/?ni  +  (-rTW  -l)sinn*;2)l  (9) 

with  the  definition 


In  general,  one  can  let  Thus,  formula  (9)  can  be  written  as 

(X 

If  an  appropriate  transport  system  is  chosen  that  can  ensure  that  the  electron  beam  forms 
a  waist  at  the  oscillator  entrance,  and  if  the  beam  radius  Rq=R^,  then,  in  later  oscillator 
transports,  the  electron  beam  will  retain  an  unchanging  beam  radius  of 

RH2)  =  RI  (12) 

As  formula  (11)  demonstrates,  the  greater  the  deviation  value  of  the  incident  beam  radius 
R^,  (when  a  is  either  much  larger  or  smaller  than  1),  the  larger  the  range  in  which  the 
electron  beam  radius  vibrates  (aR  ~  )  will  be. 

^  Ot  ^  Ot  ' 

OL 

We  calculated  the  effects  of  different  emittances  on  optical  power  where  R^  =0  and 
(see  Figure  3).  When  emittance  is  doubled,  optical  power  decreases  by  more  than 
twenty  times. 
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power  P/W 


Fig.  2  Ouput  power  at  different  normalized  electron  beam  emittances 
4.  Nonideal  incidence  of  electron  beams 

Nonideal  incidence  is  encountered  most  frequently  in  FEL  nonlinear  rotatioiL  This 
has  practical  significance  for  the  fields  of  physics  and  engineering.  One  situation  of  nonideal 
incidence  is  when  electron  beam  injection  deviates  from  the  oscillator  axis,  and  another  is 
when  the  position  of  injection  is  correct,  but  there  is  a  deviation  in  the  direction  of  injection. 
Actually,  nonideal  incidence  may  be  a  combination  of  both. 

4.1  Electron  beam  ujection  aiming  position  deviation  (Deviation  incidence) 

Before  the  electron  beam  enters  the  oscillator,  it  must  be  guided  through  the  beam 
transport  system  to  an  appropriate  position  and  have  a  definite  phase  geometry  (see  the 
previous  section).  Because  of  transport  system  mechanical  processing  errors, 
electromagnetic  wave  motion  and  outside  interference  can  bring  about  aiming  position 
deviation.  Physics  and  design  cannot  demand  that  the  above-mentioned  factor  be  zero,  but 
can  only  propose  a  maximum  allowable  deviation.  In  addition,  it  is  impossible  for  two- 
dimensional  calculation  of  rotational  symmetry  to  realize  this.  It  can  be  predicted  that 
electron  beam  deviation  incidence  will  lead  to  nonrotational  symmetiy  of  light  field 
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transverse  distribution  (this  can  be  an  aspect  of  the  light  guide  phenomenon)  and  will  have 
an  influence  on  optical  power. 
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Fig.  3  Output  power  at  different  levels  of  electron  beam  incident  deviation 

We  calculated  light  field  output  power  at  different  levels  of  electron  beam  incident 
deviation^  (see  Figure  3).  One  can  see  that  output  power  decreases  exponentially  as  the 
amount  of  deviation  increases.  When  deviation  is  an  electron  beam  radius  magnitude,  light 
power  decreases  by  two  magnitudes.  Thus,  controlling  deviation  so  that  it  is  less  than  1/10 
of  the  beam  radius,  or  about  0.3  mm  (under  the  parameters  we  use)  is  a  problem  to  which 
engineering  designers,  especially  electron  beam  transport  system  designers,  must  pay 
attention. 

Figure  4  describes  the  evolution  of  light  field  transverse  distribution  when  deviation 
is  0.2  cm.  In  order  to  better  understand  the  essence  of  the  evolution  diagram,  an  electron 
beam  envelope  evolution  diagram  with  the  same  conditions  is  provided  (see  Figure  5). 
[Under]  matching  incidence  conditions,  Rq^R^c  is  obtained.  Despite  [or]  because  of 
deviation  incidence,  the  beam  envelope  oscillates  and  wiggles  above  and  below  the  principal 
axis,  but  the  envelope  radius  does  not  fundamentally  change.  Combining  Figures  4  and  5 
[shows  that]  beginning  at  the  Im  section,  the  light  field  itself  is  smaller.  This  is  the  linear 
gain  section.  Although  the  electron  beam  deviates  [at]  the  y>0  area,  the  light  field 
maintains  left-right  symmetry  (axial  symmetry).  At  2m,  a  small  error  appears  in  the  light 
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Piancha  (deviation)  is  translated  by  the  Chinese  authors  as  "displacement." 
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field,  indicating  that  a  higher  level  mode  appears,  even  though  the  electron  beam  is  already 
oscillating.  Later,  as  the  electron  beam  wiggles  more  toward  the  y <0  area,  the  light  field 
also  deviates  towardy  <0  because  of  exponential  gain.  Finally,  a  clear  left-right  asymmetry 
is  formed.  This  shows  that  the  light  guide  phenomenon  is  primarily  located  in  the 
exponential  gain  area. 


o 


\ 


o 


o 


kj 


Fig.  4  Evolution  of  light  field  transverse  distribution  where  electron  beam  incident 

deviation  is  0.2  cm 


a.  z=0;  b.  z  =  lm;  c.  z=2m;  d.  z=3m;  e.  z=4m;  f.  z=5m 
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Figure  6.  Output  power  at  different  electron  beam  incident  angles  (in  units  of  jJ^m) 
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Figure  8.  Electron  beam  locus  where  incident  angle  is  and  i?o=^oe 
5.  Conclusions 

This  work  is  some  of  the  results  of  our  nonaxial  symmetry  2D  FEL  calculations. 
There  is  still  much  work  to  be  carried  out  in  this  area.  But  in  light  of  the  results  that  have 
already  been  obtained,  this  kind  of  nonaxial  symmetry  2D  approximation  model  not  only  has 
value  for  physics,  but  also  has  great  significance  for  engineering  design.  At  least  we  have 
suggested  some  requirements  for  pre-oscillator  electron  beam  transport  systems.  Judging 
from  the  calculated  results,  the  calculation  method  and  program  organization  of  this 
[computer]  program  are  reasonable  and  reliable.  This  program  can  provide  all  of  the 
information  about  light  fields  and  electron  beam  evolution,  (including  light  fields’  transverse 
amplitude,  phase  distribution  electron  beam  energy  spectra,  and  phase  diagrams),  and  thus 
is  beneficial  for  physics  research.  Its  greatest  difference  from  ordinary  programsf*'^!  is  that 
it  can  process  nonaxial  symmetry  effects,  such  as  nonideal  injection  of  electron  beams.  Of 
course,  this  program  can  also  be  used  for  variable-parameter  oscillators. 

Up  to  now,  we  have  only  made  some  amplifier  computations,  and  oscillator 
operations  required  approximately  190  CPU  hours  (using  an  AST  386/25  microcomputer). 
Thus,  further  reduction  work  remains  to  be  done. 
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